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Abstract 



We make progress towards proving the strong Eshelby's conjecture in three 
dimensions. We prove that if for a single nonzero uniform loading the strain in- 
side inclusion is constant and further the eigenvalues of this strain are either all 
the same or all distinct, then the inclusion must be of ellipsoidal shape. As a con- 
sequence, we show that for two linearly independent loadings the strains inside 
the inclusions are uniform, then the inclusion must be of ellipsoidal shape. We 
then use this result to address a problem of determining the shape of an inclu- 
sion when the elastic moment tensor (elastic polarizability tensor) is extremal. 
We show that the shape of inclusions, for which the lower Hashin-Shtrikman 
bound either on the bulk part or on the shear part of the elastic moment tensor 
is attained, is an ellipse in two dimensions and an ellipsoid in three dimensions. 

1 Introduction and statements of results 

In theory of composites or micro-structures, it is important to find inclusion shapes 
which produce the minimal energy. In relation to such shapes Eshelby [13] showed 
that if the inclusion is of ellipsoidal shape, then for any uniform loading the strain 
inside O is uniform. We call this remarkable property Eshelby's uniformity property. 
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Eshelby then conjectured in [T3] that elhpsoids are the only shape (structure) with 
such a uniformity property. 

Eshelby's conjecture may be interpreted in two different ways: 

Weak Eshelby's conjecture. If the strain is constant inside $7 for all loadings, 
then 17 is an ellipse (2D) or an ellipsoid (3D). 

Strong Eshelby's conjecture. If the strain is constant inside for a single 
loading, then 0, is an ellipse (2D) or an ellipsoid (3D). 

The strong Eshelby conjecture of course implies the weak one. 

The strong Eshelby's conjecture has been proved to be true in two dimensions by 
Sendeckyj [28] (see also \19\ I22j for alternative proofs). However it is only recently 
that the weak Eshelby conjecture was proved to be true in three dimensions: by 
Kang-Milton [191 ^22j. We refer to above mentioned papers (and [Lhi) for 

comprehensive account of developments on the Eshelby conjecture. 

Regarding the strong Eshelby conjecture in three dimensions, important pro- 
gress has been made by Liu. He showed in [22] that the conductivity version of 
the strong Eshelby conjecture fails to be true completely (see \22} for a precise 
statements.) However, the strong Eshelby's conjecture (for elasticity) has not been 
proved or disproved. In this paper we consider the strong Eshelby conjecture. Even 
though we are not able to resolve the conjecture completely, we obtain results which 
is stronger than the weak version of Eshelby's conjecture (and weaker than the 
strong version). We show that if the strain inside inclusion is constant and in 
addition the eigenvalues of the constant strain are either all the same or all distinct, 
then the inclusion is of ellipsoidal shape. We then use this result to show that 
for two linearly independent loadings the strains inside the inclusions are uniform, 
then the inclusion must be of ellipsoidal shape. It is worth emphasizing that the 
weak Eshelby's conjecture requires 6 linearly independent loadings while the strong 
Eshelby's conjecture does a single loadings. 

In order to present results in more precise way let us introduce some notation. 
Let be a bounded domain with a Lipschitz boundary in M*^, d = 2,3. The domain 
0, is occupied by a homogeneous isotropic elastic material whose Lame parameters 
are A and jl. We assume that the background (the matrix) is also homogeneous and 
isotropic, and its Lame parameters are A and /U. Then the elasticity tensors for the 
matrix and the inclusion can be written respectively as 

C° := AI I + 2/xI and := AI (g) I + 2/iI, (1.1) 

where I is the d x d identity matrix (2-tensor) and I is the identity 4-tensor. The 
the elasticity tensor for M'^ in the presence of the inclusion is then given by 

Cf, := (1 - ln)CO + loC\ (1.2) 

where In is the indicator function of 0. 
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Let K and k be bulk moduli of R"^ \ 17 and 17, respectively, namely, 
K = dX + 2fi and k = d\ + 2fi, d = 2,3. 
It is always assumed that the strong convexity condition holds, i.e., 

fi> 0, Hi> 0, fl> and k > . (1.3) 

We also assume that 

(A-A)(/i-^) >0, 

which implies that — C'^ is either positive or negative definite as an operator on 
the space MJ of all d x d symmetric matrices. 

We consider the following problem of the Lame system of linear elasticity: For 
a given non-zero symmetric d x d matrix A 

V • Cn£(u) =0 in M'^ 

(1-4) 

u(x) — Ax = Odxl"*" ) as |x| — > oo, 
where ^(u) is the strain tensor, i.e., 

£{u) := 2 (^'^ + Vu"^) for transpose). 

The matrix A represents a uniform loading at infinity. 

In this paper we prove the following improvements of the weak Eshelby conjecture 
for the three dimensional elasticity. 

Theorem 1.1 Let Q. he a simply connected hounded domain in R? with a Lipschitz 
boundary. If the strain tensor £{u) of the solution u to is constant in Q for 

a nonzero symmetric matrix A and £{u) within Q has eigenvalues which are either 
all distinct or all the same, then Q is an ellipsoid. 

Theorem 1.2 Let Q, be a simply connected hounded domain in M'^ with a Lipschitz 
boundary. If the strain tensors of solutions to ^.4\ ) for two linearly independent A 's 
are constant in Q, then H is an ellipsoid. 

The second main result of this paper is on the shape of the inclusion whose 
elastic moment tensor (elastic polarizability tensor) has an extremal property. In 
order to explain the second result, we take the following definition of the Elastic 
Moment Tensor (henceforth denoted as the EMT) [5, Lemma 10.3]: Let A be a 
dx d matrix and let ua be the solution to (11. 4p corresponding to A. Then the EMT 
M associated with the inclusion $7 and the elasticity tensors C*^ and C"*^ is a 4-tensor 
defined by 

MA = /" (C^ - C°)f (ua) cfx. (1.5) 
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The EMT may be defined in many different but equivalent ways. It is worth noticing 
that if the strain (?(ua) = B is constant in 0, then 



MA = IJ1|(C^ -C°)B, (1.6) 

where \Q\ denotes the volume of Q. 

The EMT enjoys several important properties. For example, it is symmetric and 
positive-definite or negative-definite on the space MJ of d x d symmetric matrices, 
depending on the sign of jl — fi. The notion of EMT is being used in variety of 
contexts such as detection of small elastic inclusions for non-destructive evaluation 
and medical imaging [H [U [3l [71 [HI [T7] and effective medium theory [5l [H [23] . 

Let us introduce more notation in order to recall the optimal trace bounds (the 
Hashin-Shtrikman bounds) for the EMT. Let 

Ai:=il(8)I, A2:=I-Ai. 
d 

Then the elasticity tensor C*^ may be written as 

C° = c/kAi + 2/xA2, 

and likewise for C"^. Since for any d x d symmetric matrix A, I 1(A) = tr (A)I 
and 1(A) = A, one can immediately see that 

AiAi = Ai, A2A2 = A2, A1A2 = 0. (L7) 

We are now able to recall the optimal trace bounds for the EMT. For d = 2,3, 



let 



K - ^ d^ + 2(d-l)/i 
^ ■ d{K- k) dK + 2{d-l)fi' ^ ■ 



d^ + d-2 ,fd-l d-1 



2fi dK + 2{d-l)fi 



(1-9) 



The following trace bounds were obtained by Lipton [2T] (see also [9]): Suppose 
|r2| = 1 and let M be the EMT associated with 0, then we have 

tr (AiM^^Ai) < Ki, (LIO) 
tr (A2M-IA2) < K2, (Lll) 

provided that k — k > 0. (If /t — k < 0, the inequalities change the direction.) Since 
AiM~-^Ai and A2M~^A2 are block diagonal components for M~^, one can see that 

trM"^ = tr (AiM^^Ai) + tr(A2M"^A2), 
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and hence 



(1.12) 



Note that AiMAi and A2MA2 are the bulk and shear parts of M, respectively. 
We also note that (jl.lOp and (jl.lip are lower bounds for M since they are upper 
bounds for M~^. It is worth emphasizing that upper bounds for M are also derived 
in [21]. In [9], it is shown that inclusions Vl whose trace is close to the upper 
bound must be infinitely thin. The upper and lower bounds for the EMT may also 
be derived as a low volume fraction limit of the Hashin-Shtrikman bounds for the 
effective moduli of the two phase composites, which was obtained by Zhikov [29^ l30] 
and Milton-Kohn [23]. Benveniste [8j obtained the upper and lower bounds of EMTs 
when those EMTs happen to be isotropic. (See also [25].) 

In this paper we are interested in the shape of the inclusion whose EMT satisfies 
the equality in either (jl.lOp or (jl.lip . This is an isoperimetric inequality for the 
EMT. In this direction we prove the following theorem: 

Theorem 1.3 Let Q be a simply connected bounded domain in M'^ with a Lipschitz 
boundary. Suppose = 1 and let M be the EMT associated with 0. // the equality 
holds in either 1^1. 10\) or 1^1. then $7 is an ellipse in two dimensions and an 
ellipsoid in three dimensions. 

We remark that optimal shapes for a cavity (hole) in two dimension were investigated 
by Cherkaev et al [10] and Milton et al [26j . 

The dimension of the space of symmetric 4-tensors in the three dimensional 
space is 21, and hence the equalities (jl.lOp and (jl.lip are satisfied on a 19 (21 — 2) 
dimensional surface in tensor space. However ellipsoid geometries (with unit volume) 
only cover a 5 dimensional manifold within that 19 dimensional space. 

It is interesting to notice similarity of Theorem 1 1.31 to the Polya-Szego conjecture, 
which asserts that the inclusion whose polarization tensor has the smallest trace is 
a disk or a ball. The Polya-Szego conjecture was proved to be true by Kang- Milton 
[18\ I19|. As for the Polya-Szego conjecture, Theorem II. 3j which concerns elasticity, 
will be proved using Eshelby's conjecture. 

In order to prove Theorem 11.31 we will show that if equality holds in (jl.lOp , 
then the strain tensor corresponding to a certain uniform loading A (with a special 
structure) is constant in 0, while if the equality holds in (jl.lip . then the strain 
tensors corresponding to five (two in 2D) linearly independent uniform loadings are 
constant in Q . Thus in two dimensions the strong Eshelby conjecture immediately 
implies that the inclusion is an ellipse. However, in three dimension, the weak 
Eshelby conjecture does not guarantee that the inclusion is an ellipsoid. In order to 
apply the weak Eshelby conjecture, we need to have equalities in both (jl.lOp and 
(jl.lip . or the equality in the whole lower trace bound (jl.l2p . But we are able to 
show additionally that if the equality holds in (jl.lOp then the eigenvalues of the 
strain tensor are all the same, and that if the equality holds in (jl.lip then strains 
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corresponding to five linearly independent loadings are constant. Thus thanks to 
Theorem 11.11 and Theorem 11.21 we are able to conclude that the inclusion is of 
ellipsoidal shape. 

This paper is organized as follows. In section 2, we show that the displacement 
vectors can be decomposed in a way similar to the Helmholtz decomposition. This 
is done using the single layer potential for the Lame system. Theorem 11.11 is proved 
in section 3, Theorem 11.21 in section 4, and Theorem 11.31 in section 4. The appendix 
is for the proof of Lemma 14.11 which is used to prove Theorem 11.21 



2 Single layer potential 

Let us first recall the notion of the single layer potential for the Lame operator 
>Ccu := V • C£'(u). The Kelvin matrix T = {Tij)^ of the fundamental solution to 
the Lame operator Cc in three dimensions is given by 



r.,(x):=-^f^-^^, x^O, (2.1) 
47r X An xr 



where 

ai = -(- + — 1 and 02 = - ( (2.2) 

The single layer potential of the vector valued density function f on dO, associated 
with the Lame parameters (A, fi) is defined by 

5n[f](x) := / r(x - y)f(y) da{y) , x G _ (2.3) 
J on 

Using the divergence theorem, we have 

5f,[f](x) = [ T^,da{y) - ^ / ^^(x - y) • i{y)da{y) 
4vr Jon |x - y| Att Jqq |x - yj'^ 

- / T^daiy) + / ^J^^ll^daiy) 

Jan |x - y| 47r Jgn |x - y| 



An 

-? / ^^V.((x-y).f(y))da(y) 
Air Jqq |x - y| 

ai +Q2 f f(y) , /■ (x-y)-f(y) 

Att 
ai + 02 



Att 

Since A|x| = 2jx[~^, we have 



Jan |x - y| 47r Jqq |x - y| 

/ ^^da(y) + ^VV • / |x-ylf(y)da(y). 
Jan |x - y| 47r Jqq 



5n[f](x) = -^i±^A / [x-y|f(y)d^(y) + ^VV. / |x-y|f(y)da(y) . (2.4) 
oTT Jqq Att Jqq 
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Let 

Wa[f](x):=^ / |x-y|f(y)da(y) . (2.5) 
Then, in summary, we have 

5o[f](x) = -^l±^A7^n[f](x) + 02 VV • Hnm^) ■ (2-6) 

It is worth emphasizing that A^7iQ[f] = 0, i.e., 7^n[f] is biharmonic, in Q and 
Thus (j2.6p shows that the solution to the Lame system in a bounded domain in Q 
or the exterior M'^ \ 17 can be decomposed into a part harmonic in 17 or \ 17 and 
a gradient part. 

Suppose that the solution u to (II. 4p inside 17 is given by 

u(x) = Bx + V, X G 17 (2.7) 

for some constant symmetric matrix B and a constant vector v. Then the solution 
is given by 

|A. + s„[fl(x), xeninn. p,, 

[Bx + V , X G 17 , 

where 

f = (Ci - Co)^(Bx)n . (2.9) 

Here n = (ni, n2, ns) is the unit outward normal vector field to 917. See [19. Section 
4]. Note that 

(Ci - Co)^(Bx)n = [(A - A)tr (B)I + 2(/u - /i)B]n. (2.10) 

Let us put 

B* := (A- A)tr(B)I + 2(/i-/u)B, (2.11) 
so that f in (j2.8p is given by 

f = B*n. (2.12) 

According to (j2.4[) . we have 

5f,[B*n](x) = -^l±^A7^n[B*n](x) + a^VV ■ Hn[B*n](x) . (2.13) 

One can easily see that 

7in[B*n](x) = B*7^n[n](x) = -B*Vpn(x), x G 17, (2.14) 
where is defined by 

po(x) :=^^|x-y|dy , x G (2.15) 
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Therefore we have 

5n[B*n](x) = ^l±^B*VApn(x) - aaVV • B*Vpn(x). 

For a 3 X 3 symmetric matrix B, let Ab := V • BV. We note that Ai = A, the 
usual Laplacian. We then define 

w;^(x) := ABPn(x), x G Ml ^2.16) 

In particular, we write wq = Wq- Then, one can easily see that 

wn{^) [ ^^da{y), x G (2.17) 

47r Jq |x - y| 

which is (2 times) the Newtonian potential of 0,. 

It is appropriate to recall now the proof of the weak Eshelby conjecture by Kang 
and Milton. In |19j . the matter was reduced to the statement: 'The Newtonian 
potential is quadratic in if and only if is an ellipsoid', which was proved by 
Dive [12] and Nikliborc }27j in relation to the Newtonian potential problem (see also 
|11]). This statement can be rephrased as 

wn is quadratic in 0, if and only if is an ellipsoid, (2-18) 

If we further put a := "^q"^ i then we have 

5n[B*n](x) = — \aB*Vwn{^) - VwE* {^)] ■ (2.19) 
We emphasize that a > 1. 

3 Proof of Theorem 11.11 

Suppose that the solution u to (|1.4|) is linear in O and given by (|2.8|) . Then by (|2.12|) 
we have 

5n[B*n](x) = (B - A)x + v, x G f). 
It then follows from ([239|) that 

aB*Vt(;Q(x) - Vw^'(x) = a2(B- A)x + a2V, x G 1). (3.1) 

Note that if eigenvalues of B are either all the same or all distinct, so are those 
of B*. After rotation if necessary, we may assume that B* is diagonal, say 

B* = diag[6i,62,63]. (3.2) 
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.abjWQ — Wq ) = linear in 17, j = 1,2,3, 

CJJL' n 



(i) Suppose first that all eigenvalues of B* are the same, i.e., 61 = 62 = 63 = h. 
In this case, since w^* = bwn, it follows from (j3.ip that 

6(a — l)Vii;f7 = linear in 17. 

Since a > 1, is quadratic in fi, and hence is an ellipsoid by (I2.18p . 

(ii) Suppose now that all eigenvalues of B* are distinct, i.e., bi 7^ hj \i i ^ j. In 
this case, (jS.ip yields that 

and hence 

abjwn - w^* ^ fji^) in n, j = 1, 2, 3, 

for some function fj which is independent of xj. Here and afterwards w denotes the 
equality up to a quadratic function. It then follows that 

^ /i - /2 B* ^2/1 - „x 

~ T r' ~ 1 ' (3-3 

oi - 62 oi - 62 

and 

(63 - 62)/! + (61 - fe3)/2 + (62 - 6i)/3 ~ 0. (3.4) 

Since fj is independent of Xj for j = 1,2,3, one can easily see that (13. 4|) holds 
only when fi, f2 and take the form 



/i(x) P 


^ m(j;3) - 


- n{x2) 


63- 


-62 


/2(X) ^ 


^ - 


m{x3) 


bi- 


63 


/3(X) P 


^ n{x2) - 


r(a;i) 


62- 


61 ' 



for some functions m, n and r. It then follows from (13.31) that 



_ m(x3) ^^(3:2) ^(a^i) . ^ 

^ (53 _ - 63) ^ (62 - b,)ib3 - 62) (61 - 63)(62 - 61) ■ ^ ■ ^ 

Since Awq = 2 in 17, we have 

m"{x3) n"{x2) , r"{xi) 

+ t; r^TT, = constant. 



(63 -62) (61 -63) ib2 - bi)ib3 - b2) (6l-63)(&2-6l) 

Thus r, n, and m are quadratic functions of xi, X2, and X3, respectively, and hence 
is quadratic in 17. Thus 17 is an ellipsoid. 

This completes the proof. □ 
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4 Proof of Theorem 11.21 

In order to prove Theorem ll.2l we use the following lemma whose proof will be given 
in the appendix. 

Lemma 4.1 Let Bi and B2 he two symmetric 3 x 3-matrices. If Bi + tB2 has a 

multiple eigenvalue for all real numbers t, then Bi and B2 can be diagonalized by 
the same orthogonal matrix. 

Let Ai and A2 be two linearly independent symmetric 3x3 matrices and suppose 
that solutions ui and U2 to ()1.4p with A = Ai and A = A2 are linear in 
Put Bj := £{uj), j = 1,2. Since the EMT is positive or negative definite on 
M| (Theorem 10.6 of [5j), (jl.6p shows that Bi and B2 are linearly independent. 
According to (|2.12|) we have 



5n[B^n](x) = (Bi - Ai)x + vi, x e J], 
5o[B*n](x) = (B2 - A2)x + V2, x e f]. 



It then follows from ([239|) that 

jaBlVwn{x) -Vw^'-ix) = a2(Bi - Ai)x + a2Vi, 
\aB2Vwn{x) - Vw^'^{x.) = Q2(B2 - A2)x + 02 V2, 



(4.1) 



for X G $7. 

Let us suppose that all of Bi, B2, and Bi + tB2 {t G M) have an eigenvalue of 
multiplicity 2 (otherwise we apply Theorem 11.11 to conclude that is an ellipsoid). 
By Lemma [4. H Bi and B2 can be diagonalized by a single orthogonal matrix. Thus 
we may assume that Bi and B2 are diagonal. Then from (j2.1ip B^ and B2 are also 
diagonal and we may let 

B* = diag[6i,6i,ci], = diag[62, &2, C2], 

where bi 7^ ci and h2 ^ C2- Since B\ and B2 are linearly independent, we have 

biC2 / C162. 



By ()4.ip . we have 



abiwn - w^^ ^ fixs), (4.2) 

B* 

aciwn - w^^ ^ g{xi,X2), (4.3) 

B* 

ab2Wn -w^^ ^ h{x^), (4.4) 

B* 

ac2W^ - w^^ ^ l{xi,X2), (4.5) 
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for some functions /, g, h, and /. Here again ^ denotes the equality up to a quadratic 
function. By (|2.16p . we have from (fOj) and ([0|) 



and from (lOI) and (Ii3]) 



(a - 1)62/(2:3) - (a - 1)61 /i(x3) s::: (a - l)(5iC2 - ^2Ci)-^^, (4.6) 



(62 - ac2)s'(xi,X2) - (61 - aci)/(xi,X2) « (1 - a)(6iC2 - 62Ci)-^^. (4.7) 



It then follows that 

We then obtain from ig^-^^ that 

(a - 1)61 w;f7 ~ /(a^a), (a - l)62';«n ~ /i(a:;3), 

and 

(aci - 6i)t(;Q g{xi,X2), (ac2 - 62)';«n ~ l{xi,X2). 

Thus we conclude that Wfi ?a 0, and hence J7 is an ellipsoid. 

This completes the proof. □ 



5 Proof of Theorem 11.31 

The space MJ is equipped with the inner product A : B, where A : B denotes the 
contraction of two matrices A and B, i.e., A : B = "^ijCiijbij = tr(A^B) where 

tr(A) denotes the trace of A. For d = 2,3, let d^: := ^^i^hil^ which is the dimension 
of AfJ. Let Bi = -^12 be a basis for Ai(MJ) (of a unit length), and {B2, . . . , B^^} 
be an orthonormal basis for A2(MJ). Then {Bi, . . . ,Brf^} is an orthonormal basis 
for MJ, i.e., 

Bi . Bj — 6ij , 

where 5ij is Kronecker's delta. Note that for any symmetric 4-tensor T, we have 

trT = ^TBfc:Bfc. (5.1) 

k=l 

We deal with the case when — C*^ is positive definite so that M is a symmetric 
positive-definite linear operator on MJ. The other case can be treated in the exactly 
same way. 
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Let us first invoke some facts proved in [9j. Introduce a 4-tensor M by 
A : MA = min / (^:(v) + l^GA) : (£'(v) + l^GA) dx 



(5.2) 



+ \n\ (C^ -C°) (C^) ^C°A : A 



for A G MJ, where 

G :=I- (C^)-^C°. (5.3) 



Note that the minimum in ()5.2p is attained by v = u — Ax, where u is the solution 
of dH])- It is proved in % Corollary 3.2] that 

AiMAi = AiMAi and A2MA2 = A2MA2 . (5.4) 

In particular, we have 

tr(AiM~^Ai) = tr (AiM^^Ai) and tr (A2M"^A2) = tr (A2M"^A2). (5.5) 

Let C be an isotropic 4-tensor, i.e., 

C = AI O I + 2/xI = dKAi + 2/iA2, 

for some A and /i satisfying (jl.3|) . Let (M°',MJ) be the space of square in- 
tegrable functions on valued in MJ and i7^(M'^,MJ) the Sobolev space. For 
P G L2 (R'^,MJ), we define Fc(P) by 

Fc(P) := -f^cHV • P), (5.6) 

where Cc = V • C£. In other words, if $ is a unique solution in H'^{W^,Aq) to 

/:c ($) + V • P = 0, (5.7) 

then Fc(P) is given by 

Fc(P) =^(1>). 
If $ is the solution to (15.71). then 



for ah ^ G H'^{R'^, MJ), and hence by taking = ^ we have 

/ P : Fc(P) = - / CFc(P) : Fc(P) (5.8) 

We prove Theorem 11.31 using the following two propositions whose proofs will be 
given at the end of this section. 
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Proposition 5.1 Let 

Sa(C°,P)= f P : Fco (If^P) + / (C°-Ci)"^P:P + 2 / P:A. (5.9) 

Jo, Jn Jn 



Then the following holds 

A:MA= sup Sa(C°,P). (5.10) 

PeL2(Md.7v/j) 

Furthermore, this supremum is attained by P = 1q (^C^ — C^) £{u), where u is the 
solution of 



We then show that structures reaching the lower trace bounds have a particular 
structure, as explained by the proposition below. 

Proposition 5.2 // equality in il.lO\) holds, then we have 

SA(CO,lnBi)= sup -Ea(C°,P) (5.11) 

PeL2(Md:MJ) 

with A = M^^Bi. // equality in lll.ll]) holds, then we have 

EA(C°,lnBfc)= sup ^a(C°,P) (5.12) 

PeL2(iRd.7v/|) 

with A = M~^Bfc, fork = 2,...,d^. 

Proof of Theorem \1. 31 Introduce a bilinear form J^^o{Q, R) by 

.Fco(Q,R)= / Q:Fco(lnR)+ / (C^ - C^)"' Q : R . 
Jn Jn 

It follows from (15.81) that 



-^co(laQ, laQ) = - / COFco (l^Q) : Fco (InQ) - [ (C^ - C°) ' (l^Q) : (InQ) 
Jn Jn 



<_ / (ci-CO) '(lcQ):(lnQ) 



f7 



< -K\\lnQ\\L2i 



for some positive constant K. The last holds due to the positive-definiteness of 
(j-<i _ (j-<o_ jyg consequence, J^£o is negative definite when restricted to H = {P £ 
L2(0 : MJ)2, supported in O}. Therefore Sa(C°,Q) = JPco(Q.Q) + 2 Q : A is 
a strictly concave functional on H, and therefore admits at most one maximizer in 
H. 
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We observe that since C — C is isotropic, if B is diagonal and all the eigenvalues 
are the same, so is (C^ — C'^)~^B. If B is trace- free and all the eigenvalues are distinct, 
so is (C^ -C°)-iB. 

Suppose that equality holds in (jl.lOp . It then follows from Proposition [5TT1 (j5.2p . 
and uniqueness of the maximizer in 17 that 

-C°)£{ui) = Bi in n, 

where ui is the solution to ()1.4p with A = M^^Bi. Recall that Bi = -^I- Therefore, 
£^(ui) is constant in Q and all the eigenvalues of <f (ui) are the same. Thus Q is an 
ellipse or an ellipsoid due to Theorem II. li 

Suppose now that equality holds in (|1.11[) . Then for similar reasons we can 
deduce that for each k = 2, . . . , d^, 

(C^ - C°) £{uk) = Bfe in n, 

where is the solution to (II. 4p with A = M^^B^. Thus is an ellipse or an 
ellipsoid due to Theorem 11.21 

This completes the proof. □ 

Proof of Proposition I5.il Following the notation of [9], we define WAiC,P), for 
A G MJ, by 

Wa{C,P)= [ P:FcP+[ (C-Cn)"^P : P+2 /p : (C^ -C)"^ (C^ -C°) A. 

It is proved in [9, Proposition 4.1], following the variational strategy given in |20) 
for the derivation of Hashin-Shtrikman type bounds, that for any isotropic elasticity 
tensor C < C''(< C^) we have 

A:MA = A: (C^-C°) (C-C^)~^(C-C°) A+ sup Wa(C,P). (5.13) 

PeL2(R'i:Mj) 

Note that the supremum is attained by 

P = ln - C") A + {Cn - C) (^(u) - A) , (5.14) 

where u is the solution to ()1.4p with A in above identity. Since (C^ — C^)(C — 
C^)~^(C — C'^) is positive definite, by sending C to C*^, and restricting the supremum 
to fields P such that l^P = P we obtain 

A : MA > sup £;a(C°,P) . (5.15) 

P6L2(IRd,AfJ) 

For any P E (M'^,MJ), and any positive definite isotropic elasticity tensor 
C < CO, we define Sa(C,P), for A € MJ, by 

Sa(C,P)= I P:Fc(lnP)+/" (C - Ci)"^P : P+2 [ P : (C^ - C)"^(C^ - C°)A. 
Jn Jn Jn 
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Note that this definition is consistent of that of Ea{C^, P) given in (j5.9p by passing 
to the limit in C. 

Introducing the decomposition P = Pq + Pu, with Pfjlfj = Pq, and Pu^n = 0, 
we have 

TyA(C,P) = EA{C,Pn) + WA{C,Pu) 

+ [ Pn : FcPu+ [ Pu ■■ FcPn- 

Let C = Co — el, where e > 0. Then we have 

Wa{C,P) = Sa(C,Pc) - e"'||Pf/|li2(Kd) + R('Pu,Pn), 

where 

R{Pu,Pn):= [ Pn:FcPu+[ Pu-FcPn+[ Pu ■ FcPu- 

jRd JiRd J^d 

By integration by parts, and by the Cauchy-Schwartz inequahty, we readily obtain 
that for e small enough R{Pu,Pfi) satisfies 

\R{Pu,Pn)\ < i^llPt/ll^a^Kd) (||Pt/||i2(-Kd) + ||Pn|li2(]Kd)) , 

where the constant K is independent of Pu, Pq, and e. As a consequence, for e 
small enough, 

-e-^\\Pu\\l,^^,-^+ R(Pu,Pn) <SKe\\Pn\\l,^^,y 

Note that from (jS.Sp JqPq '■ Fc (Pn) is negative definite, therefore 

EAiC,P) < K||Pn|li2(]R<.)(-||Pn||i2(K.) + 1), 

where K is another constant independent of Pn and e. Thus ||Pn||/,2(^iRd~) must stay 
bounded, uniformly with respect to e, close to the supremum. Taking the limit as e 
tends to zero we obtain (j5.10p . Replacing C by C'^ in (j5.14p concludes the proof. □ 

Proof of Proposition \5.'A Given k G {1, . . . , d^:}, choose A = M^^Ai(Bfc) and use a 
test function P = l^Ai(Bk) in (|5.10p . This gives 

M-^AKB,,) : Az(Bfc) > t^A(C°, InM^k)) (5.16) 
= / A;(Bfc) : Fco (InA^B,,)) + / (C" - C^)-' A^B,) : MB^) 

+ 2 / A;(Bfc) :M-iA,(Bfc). 
Jn 
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Summing these inequalities over k, we obtain 

tr (A^M-^AO > V / Ai{Bk) : Fc« {InM^k)) (5.17) 

+ E / (C°-Ci)"'Ai(Bfc) :AKBfc) + 2tr(A;M-%). 



It is proved in [9, (4.27) & (4.28)] that 



k=l 

and 



5:/^A,(B,):Fco(l.A,(B,)) = --^^ 
V / A2(Bfe):Fco(lcA2(B,))= ^ " ^ " ^ 



fc=i 

Since 

1 . 1 



(i(A + 2/i) 2/x 



d(K — k) 2(/x — /i) 

one can immediately see that 

^l^-lA,(B,):Al(Bfc^ ^ 



V / (C° - C 



(i(K — k) 



and 



V / (C°-Ci) 'A2(Bfc):A2(Bfc) 



2(/x - /i) 



fc=i ■ 

Therefore, we get 

V [ [ Ai(Bk) : Fco (InAKBfc)) + / (C° - C^)"' A^B^) : A^B^)! = -Ki 

for / = 1,2, where Ki is given in (|1.8p and (|1.9p . It then follows from (|5.17|) that 

tr (A/M"^AO > -Ki + 2tr (A^M^^A;). (5.18) 

Suppose that equality in (jl.lOp holds. Then, in view of (j5.18p . the inequality in 
(j5.17p becomes an equality, and so does the one in ()5.16p . Since Ai(Bi) = Bi and 
Ai (Bfc) = for A; = 2, . . . , , we have 

^A(C°,lnBi)= sup Ea{C^,P), A = M-1Bi. (5.19) 

PeL2(Rd:Mj) 
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Likewise, if equality in (jl.lip holds, then 

^A(C°,lf7Bfe) = sup £;a(C°,P), A = M-iBfc, (5.20) 

PeL2(R<i:A/|) 

for A; = 2, . . . , (i^K, and the proof is complete. □ 

A Proof of Lemma 14.11 

Proof. By considering jBi + B2 and taking the limit t — > 00, one can see that B2 
also has a multiple eigenvalue. If B2 has an eigenvalue of multiplicity 3, then B2 is 
a constant multiple of the identity matrix and hence the conclusion of the lemma 
holds trivially. 

Let us assume that B2 has an eigenvalue of multiplicity 2. Note that, for any 
real number s and orthogonal matrix U, 

UBiU"^ + tUB2U"^ - si = U(Bi + tBa - sI)U"^ 

has a multiple eigenvalue regardless of t. Therefore we may assume that B2 takes 
the form 

/O 0^ 
B2 = 

VO 1; 



while Bi is arbitrary, say 



Bi 




Let T{t) be the discriminant of the characteristic polynomial of Bi + tB2. Since 
Bi + iB2 has a multiple eigenvalue for all t, T(t) = 0. Then a straightforward 
calculation shows that the coefficient of term of T{t) is given 
and hence we have 

a-b = d = 0. 

It then follows that the coefficient of t"^ term of T(t) is given by (e^ + /^)^, and hence 

e = / = 0. 

Thus Bi takes the form 

Bi = 




This completes the proof. □ 
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